New universality of metal-insulator transition in integer quantum Hall effect system 
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A new universality of metal-insulator transition in integer quantum Hall effect (IQHE) system 
is studied based on a lattice model, where the IQHE states only exist within a finite range of 
Fermi energy in the presence of disorders. A two-parameter scaling law is found at the high-energy 
boundary where direct transitions from high IQHE states to insulator occur. We find p xx = p xy at 
the critical point whose value can continuously vary as a single function of the Landau-level filling 
number n„. Such a new universality well explains recent experiment by Song et al. (Phys. Rev. 
Lett. 78, 2200 (1997)). 
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The global phase diagram proposed by Kivelson, Lee 
and Zhang @ (KLZ) for the quantum Hall effect (QBE) 
system has stimulated a series of recent experiments. Al- 
though many of these experiments arc in support of the 
overall picture of KLZ, controversies arise concerning the 
select rule for metal-insulator transitions in such a sys- 
tem [0-0. The global phase diagram predicts that the 
transition from an IQHE state to insulator can only oc- 
cur at low filling number corresponding to the v = 1 
IQHE state to insulator (1 — > 0) transition, while at 
higher filling number the only allowed transitions are the 
nearest-neighbor plateau-plateau transitions (y — > v±\). 
Experimentally, however, direct transitions from v = 1, 
v = 2, v = 3, and v = 6 IQHE states to insulator have all 
been observed recently HQ] on the high-filling number 
side. Although the apparent inconsistency of the v = 2 
IQHE state to insulator transition with the global phase 
diagram may be reconciled by assuming that the low- 
est Landau level (LL) is spin degenerate ||, the higher 
plateaus to insulator transitions (3 — > and 6 — > 0) can 
not be explained similarly due to the fact that these tran- 
sitions must take place at higher LL's. 

The select rule of the global phase diagram is based 
on the so-called float-up picture where extended 

states in the weak field limit are assumed to shift to- 
wards higher energy indefinitely without disruption. But 
recent numerical calculations have shown that the 
float-up picture is not correct in a lattice model in which 
extended states merge and disappear before they can 
reach to the band center. This non-float-up picture holds 
down to very weak magnetic field limit 11 which means 
that the lattice model is not equivalent to the continuum 
model even in the limit where the magnetic length scale 
is much larger than the lattice constant. This surprising 
result, which can be well understood based on a topolog- 
ical Chern number description |J, implies that the lat- 
tice effect will always remain to be an important factor 
in the problem of metal-insulator transitions even in the 
weak field limit. It may thus provide an explanation for 
the aforementioned experimental observations of a direct 
transition from the higher IQBE states to insulator. 

Recent experimental measurement has further indi- 



cated [Q that the direct transition from a higher IQHE 
state to insulator may belong to a new universality. In 
the corresponding critical regime, the longitudinal and 
transverse conductances generally satisfy a xxc — a xyc 
within the experimental error and their value varies con- 
tinuously with the filling number, which is different from 
the usual plateau-plateau transitions (including 1 — > 
transition as well) in which a xxc — 0.5e 2 /h and a xyc = 
(y + 0.5)e 2 /h [[To). So far such a new class of metal- 
insulator transitions in the IQHE system cannot be un- 
derstood by any existing theory of the IQHE based on 
the continuum model. On the other hand, since the di- 
rect transitions from higher IQHE states to insulator are 
present naturally in a lattice model as mentioned above, 
it would be very desirable to examine the corresponding 
critical behaviors. 

In this paper, we present a numerical study based on 
the calculation of longitudinal and transverse conduc- 
tances and a finite-size scaling analysis near these critical 
points. Our results show that a direct metal-insulator 
transition in this regime indeed belongs to a new univer- 
sality. In particular, the scaling behavior of localization 
length near each critical point obeys a real two-parameter 
scaling law, in contrast to the plateau-plateau transi- 
tions whose critical behavior can be reduced to a one- 
parameter scaling law JTT[ | . As the consequence, a xxc and 
o~xyc, which are equal to each other, continuously vary as 
a function of a single parameter — the LL filling factor 
n„, and does not depend on magnetic field or disorder 
strength independently. Since such a critical behavior 
is independent of the strength of magnetic field, we are 
able to make a close comparison of our numerical results 
with experimental measurement [0 and find a consistent 
explanation. 

The lattice Hamiltonian to be studied is a tight-binding 
model (TBM) given as follows: 



H 



<ij> 
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where the hopping integral is taken as the unit, and cf 
is a fermionic creation operator with < ij > referring to 
two nearest neighboring sites. A uniform magnetic flux 



per plaquette is given as 4> = J2a a ij — 2n/M, where 
the summation runs over four links around a plaquette. 
u>i is a random potential with strength \wi\ < W/2, and 
the white noise limit is considered with no correlations 
among different sites for u>i. 

In the weak disorder limit, a well-defined IQHE plateau 
structure is exhibited in the Hall conductance. For our 
purpose, the disorder strength W is continuously in- 
creased to identify a direct transition from the ^th IQHE 
state to insulator. For instance, in Fig. 1(a), by fixing 
the Fermi energy near the center of the v — 2 plateau 
(Ej — —2.75) at weak disorder, we show how the Hall 
(cr xy ) and longitudinal (a xx ) conductances evolve with 
disorders (at a flux strength <f> — 27r/16). Here a xy is 
calculated based on the Kubo formula and a xx by the 
Landauer formula jL2j. At weak disorder (W < 1), a xy 
is well-quantized at ve 2 jh with v = 2. With W in- 
creasing from 1 to 8, a xy continuously reduces to zero 
without showing a v = 1 plateau, indicating a direct 
v = 2 — > v = transition. Two curves for a xy at dif- 
ferent sample sizes (N = 16 2 and N — 32 2 ) cross at 
a critical point which also coincides with the maximum 
point of the longitudinal conductance <r xx . The latter is 
also sample-size independent at the corresponding criti- 
cal disorder strength W c = 3.5. But away from W c , a xx 
monotonically decrease with the increase of sample size 
as a typical behavior in localized states. Thus the exis- 
tence of a single W c identified from both a xy and a xx con- 
firms that the v = 2 — > transition is really a one-step 
transition. Furthermore, a calculation of the finite-size 
localization length by using the transfer matrix method 
p3| also verifies the same W c : at the fixed Ef — —2.75, 
we find that the localization length satisfying the scaling 
relation, Xl/L = f'(£(W)/L), for a stripe sample with 
the width varying from L = 16, 24, 32, 48, 64, to 80. Here 
£ is identified as the thermodynamic localization length 
which diverges at W c = 3.5 as shown in Fig. 1(b). At 
such a critical point, we obtain a xxc = (1.06 ± 0.02)e 2 //i 
and a xyc = (1.01 ± 0.06)e 2 //i, or, cr xxc = cr xyc within the 
statistics error bar (6%). 

Fixing the Fermi energy Ef near the center of the vth 
plateaus at weak W, we always find a similar direct tran- 
sition to insulator [y — > 0) as W increases. The critical 
disorder W cv for v — > transition is found to satisfy a 
sequence W cv < W cv i for v > v 1 . For instance, we have 
W cl = 3.8, W c2 = 3.5, Wc3 = 3.2 and W c4 = 2.8 at a flux 
strength 27r/16. Up to the weakest flux strength 27r/384 
that we can get access to, we always find the same se- 
quence for the disappearance of the IQHE plateaus. Ac- 
cording to this sequence, near the critical disorder W cv 
for i^th plateau, all higher plateaus have already disap- 
peared while lower ones are still robust. We can then 
define the transition of v — » as the boundary of the 
IQHE regime on the high-energy (or high-filling number) 
side, beyond which the IQHE states no longer exist up 
to the band center. The corresponding critical conduc- 



tances a xxc and a xyc at W cu always satisfy the relation 
&xxc — &xyc within the numerical error bar, and their 
values continuously vary with the plateau index v and 
the Fermi energy Ef. Nevertheless, we find that all of 
them can be scaled into a single curve as the function of 
the LL filling number n v shown in Fig. 2, in which the 
magnetic flux changes from 27r/16 to 27r/96. In Fig. 2, 
we show p xxc and p xyc instead of a xxc and a xyc in order 
to compare with experimental data of Song, et al. j| in 
the insert of Fig. 2. We see that the overall agreement 
between the theory and experiment is quite reasonable. 
(Note that in the insert the top horizontal axis is carrier 
density which has not been converted to n v here as the 
corresponding magnetic fields are not provided in Ref. 
J!.) Both p xxc and p xyc monotolically increase with the 
decrease of n v and are eventually saturated at h/e 2 . We 
can make a closer comparison: in Ref. 0] the v = 3 — * 
transition covers a range of critical resistivities from 0.27 
to 0.4 in unit of h/e 2 , which agrees well with the theo- 
retical values shown in Fig. 2. The theory also predicts 
that the transitions reported in the experiment with crit- 
ical resistivities less than 0.27 (h/e 2 ) actually come from 
higher (y > 4) IQHE states to insulator transitions. 

Thus we have identified a new universality of metal- 
insulator transition at the boundary between the IQHE 
regime and insulator on the high filling-number (energy) 
side in the TBM. It is quite different from the well-known 
plateau-plateau transition (y — > v ± 1) within the IQHE 
regime as well as the transition from the IQHE state to 
insulator on the low filling- number side (i.e., v = 1 — > 0) 
where, at v — * v + 1 transition, one has a xxc ~ 0.5e 2 //i, 
while a xyc = (0.5 + v)e 2 jh |To| ] . It provides a unique 
explanation for recent experiments [^|^] . In order to fur- 
ther understand this new universality of critical behavior 
at the high-filling boundary of the IQHE, we investigate 
the scaling behavior in this regime below. 

According to the general scaling theory of the QHE 
system [Q, the finite-size localization length may be 
written as a general function of two parameters, i.e., 
Xl/L = f(L/£,p), at a large sample with a width L. 
Here the first parameter is L divided by the thermody- 
namic localization length £ as usual |r| , and the second 
one, p, may be chosen as either a xx or a xy at a fixed 
sample size Lq. In the present case, it is convenient to 
choose <r xxc as p (which can be connected with a xx and 
a xy at L = Lq by the scaling- flow diagram Jl4[ ). Based 
on this scaling hypothesis, the numerical data at different 
flux strength 2ir/M's, disorder strength Ws, sample size 
L's, etc., should all collapse into a single curve so long 
as the critical conductance a xxc is the same. Indeed, 
as shown in Fig. 3, the data of finite-size localization 
lengths can be well fit into one curve by choosing a sin- 
gle scaling variable L/£ at the given a xxc = 1.06. Such a 
(T IIC corresponds to a v — 2 — > transition. By changing 
the plateau index v or the Fermi energy Ef, the critical 
conductance a xxc can change continuously as discussed 
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before. One expects the scaling curve to change corre- 
spondingly. In the insert of Fig. 3, scaling curves at two 
critical points {y = 2 — ► and is = 4 ^ 0) are shown 
with a XX 2 = 1-06 and a xx ^ = 2.1, respectively. Different 
scaling curves have also been obtained for the transitions 
of v = 3 — > and v = 1 — » in this regime. All of these 
different cases can be generally specified by a single LL 
filling number n u , since a xxc is uniquely determined by n v 
as shown in Fig. 2. Therefore, a two-parameter scaling 
law is well established here which characterizes a new uni- 
versality of metal- insulator transitions between the IQHE 
regime and insulator at high-filling-number boundary. It 
is noted that the transition from the IQHE state to insu- 
lator at low-filling-number boundary (corresponding to 
v = 1 — > 0) is usually well described by a one-parameter 
scaling law |ll| , which may be understood as the second 
parameter a xxc always remains a constant in this regime 
|Tof . So the present case is the first one in the QHE sys- 
tems that a two-parameter scaling law becomes necessary 
in order to describe the critical behavior. Finally, if one 
plots the finite-size conductance a XX L as a function of 
the finite-size localization length Xl/L at different flux 
strengths (27r/16 — 27r/96), as well as different plateau 
index z/s, the Fermi energy Ef's and disorder strength 
Ws, all the data also collapse into one curve as shown in 
Fig. 4. (The notations of the data in Fig. 4 are similar 
to those in Fig. 3.) It means that u xx l is uniquely de- 
cided by Xl/L, and thus also satisfies the two-parameter 
scaling law. 

Finally, we would like to point out that in the above 
discussion, we have not considered the spin degree of free- 
dom. For the conventional plateau-plateau transition, it 
is well-known jl5| ] that the mixing of LL's for different 
spins by a spin-orbit coupling will not change the uni- 
versality class and such a spin-coupled system behaves 
still like a spinless system. In our case, we have checked 
that introducing a weak spin-orbit coupling effect will 
also make the system with spin degree of freedom be- 
have like a spinless system. For example, the v = 2 — > 
and ^ = 3^0 transitions remain to be direct transitions 
(even with Zeeman splitting) with the same universality 
of critical behavior discussed above for the spinless case. 
Both diagonal and Hall conductance are also equal to 
each other and fall into the same range shown in Fig. 2. 
Details will be presented elsewhere. 

In summary, we have identified a new two-parameter 
scaling law in the critical regime of direct transitions from 
the high IQHE states to insulator in the TBM. Such a 
new universality of metal-insulator transitions predicts 
that p xxc — Pxyc and their value depends solely on the 
LL filling number n v . The critical behavior of this lat- 
tice model in the aforementioned direct phase transition 
regime provides a consistent explanation for recent ex- 
perimental measurements. 
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Fig. 1 (a) The evolution of a xx and a xy with disorder 
strength W at a fixed Fermi energy Ef = —2.75, which 
corresponds to the center of the v = 2 IQHE plateau 
at weak W limit. Lattice sizes vary form 16 x 16 (*), 
32 x 32 (+), to 64 x 64 (•), and the sample-independent 
fixed point at W c = 3.5 indicates the crossing of the last 
extended level at Ef. (b) Thermodynamic localization 
length £ also diverges at W c . 

Fig. 2 Critical resistivities, p xxc and p xyc , as function 
of Landau- level filling number n„, respectively, at flux 
strength f§ (A, +); §f (o, x); §§ (□, *) and §£ (., 
none). The insert: the experimental data by Song, et al. 

§• 

Fig.3 
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Two-parameter Scaling function Al(L/£, a xxc )/ L versus 
L/fl at <r xxc — 1.06. All the data at different disorders 
collapse into one curve with the stripe sample width L 
varying from L = 16 to 80 and flux strengths from ^ 
(*), If (<>)> If ( D ) to |f ( A )- Thc insert: two different 
scaling curves for the v = 2 — > (at (j xx2 — 1.06) and 
v = 4 — > (at cr a;a ;4 = 2.10) transitions, respectively. 

Fig. 4. Finite-size conductance ct^l as a single func- 
tion of Xl/L at different sample-sizes, flux strengths, 
plateaus indexes, disorder strengths, and Fermi energies. 
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Fig. 2 




Fig. 3 



